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Abstract 

We study non-equilibrium quantum dynamics of the single-component scalar field theory in 1 + 1 
space-time dimensions on the basis of the Kadanoff-Baym equation including the next-to-leading- 
order (NLO) skeleton diagrams. As an extension of the non-relativistic case, we derive relativistic 
kinetic entropy at the first order in the gradient expansion of the Kadanoff-Baym equations. The 
derived entropy satisfies the H theorem. Next we perform numerical simulations in spatially ho- 
mogeneous configurations to investigate thermalization properties of the system by evaluating the 
system entropy. We find that at later times the kinetic entropy increases approaching the equilib- 
rium value, although the limited time interval in the early stage invalidates the use of it. 

1 Introduction 

Non-equilibrium quantum field theories provide a suitable framework to investigate a large variety 
of topical problems in high energy particle physics, astrophysics, cosmology, as well as condensed 
matter physicspQ [2]. In the context of heavy ion collision physics, the early time evolution of the 
colliding system toward the quark-gluon plasma (QGP) state has attracted a lot of theoretical 
interests for recent years. Success of ideal hydrodynamic models for describing bulk properties of 
the matter created at Brookhavcn's Relativistic Heavy Ion Collider (RHIC) seems suggesting that 
the produced system is strongly interacting and nearly thermalized within a short time[3] compared 
with perturbative analysis [J. There are various theoretical studies on the possibility for this short 
time thermalization, some of which rely on the instabilities in the plasma[Sl [SJ [7J |H1 GO EH HI], an< ^ 
some others include the 2-to-3 processes in parton cascade simulations [T2]. 

In the earliest stage of the high-energy nuclear collisions, the system will be so dense that it 
would be more suitable to describe the system in terms of the quantum field degrees of freedom 
than in the particle basis. As a first step of this approach toward the early time dynamics of the 
nuclear collisions, we study here the non-equilibrium Xcf) A scalar field theory in 1+1 dimensions on 
the basis of the Kadanoff-Baym (KB) equations. 

As early as 1960s, based on a functional formulation of Luttinger and Ward [13], Baym and 
Kadanoff studied the Dyson-Schwinger equation for the two-point function G(x, y) [14] . Then Baym 
reformulated it in terms of variational principle, introducing the so-called ^-derivable approxima- 
tion [T5l [T6] . The functional $[G] is given by a truncated set of closed two-particle irreducible (2PI) 
diagrams, and generates the driving terms of the equations of motion. The main virtue of this ap- 
proximation is that the resulting KB equations conserve the energy and momentum of the system. 
This approach was extended to relativistic systems and formulated using the path integral by Corn- 
wall, Jackiw and Tomboulis in [17] . It can be extended further to more general non-equilibrium 
many-body systems based on the Schwinger-Keldysh real-time path integral method [18| [T9]. 

In the last several years, the real-time field dynamics has been newly investigated by several 
authors. A seminal work was carried out by Danielewicz 20 , who for the first time studied the 



full KB equations in the context of the heavy ion collisions at non-relativistic energies. He used a 
spatially homogeneous initial condition with the non-spherical Fermi distribution for the nucleon 
momentum. Thermalization problem in the relativistic \4> 4 scalar field theory in 1 + 1 [21] [22] ,2 + 1 
[2"5] and 3 + 1 [55] dimensions, has been investigated with keeping the NLO skeleton diagrams 
in <£. Extension to the O(N) theory at the next-to-leading order in 1/N expansion can be found 
in Ref. [26] . Importantly, all these analyses indicate that thermalization is achieved in course of 
the time evolution of the system independently of the initial conditions. The number distribution 
functions of the quasi-particles are found to approach the Bose distribution. 

The approach to the equilibrium state will be quantitatively characterized if system entropy 
can be introduced properly. In fact it is an open question how to choose the gross variables and 
define the corresponding entropy of the system in general non-equilibrium situations. There is 
no entropy production in fully microscopic calculations. We use the variable G(x, y) in the KB 
approach. In the non-relativistic case, the kinetic entropy is introduced at the first order of the 
gradient expansion in Refs. M and 28 0. 

To our knowledge, the entropy production has not ever 
been estimated in the relativistic KB dynamics. Here we shall extend the entropy to the relativistic 
case in the first order of the gradient expansion. This will provide us, for example, of a criterion 
how much each microscopic process contributes to thermalization of the system. 

The KB equations deal with the evolution of the field (<fi) and the two-point function G(x,y), 
and effectively contain particle number changing processes such as 1 H 3 even in the NLO, if 
interpreted in the particle basis. In contrast, the Boltzmann equation includes only the 2 <-> 2 
scattering processes to this order, which preserves the total particle number. This difference should 
be reflected in the behavior of the system evolution, especially in the entropy production. We 
expect that this aspect of the KB equations is important to understand the possibility of the rapid 
thermalization. 

For demonstration, we shall numerically solve the non-equilibrium dynamics of \<fr theory on 
the basis of the KB equations. In order to reduce the numerical cost, we restrict our simulations 
to the spatially uniform case without the mean field (</>) = in 1 + 1 dimensions. We start the 
simulations with the non-thermal initial conditions, and show the time evolution of the system 
through the particle number distribution functions, the energy content, the entropy production 
and so on. 

This paper is organized as follows. In Sec. 2 we briefly review the formulation of the KB 
equation for the relativistic scalar field theory, using the Schwinger-Keldysh formalism and 2PI 
effective action in the NLO. Next, we present the derivation of the entropy for the relativistic KB 
equations in the first order in the gradient expansion in Sec. 3, which is the main theoretical part of 
this work. The expression for the entropy is found to be a natural extension of the non-relativistic 
one given in [271 12"51 121)1 1313] in the local approximation. This entropy satisfies the H-theorem. In Sec. 
4 we show the numerical simulations of the KB equations. The particle number distribution and 
the entropy of the system are calculated in terms of the numerical solutions of the KB equations. 
Finally Sees. 5 and 6 are devoted to discussions and summary of this study. 

2 Kadanoff Baym equation 

We briefly review the derivation of Kadanoff-Baym equation and fix our notationspQ. For the scalar 
field theory C — hd^d^cf) — ^m 2 cf> 2 — jjXip 4 , the 2PI effective action with vanishing mean field 
{4>) = (unbroken phase) is written as 

T[G] = ^Trln(G)- 1 + '-G, 1 G + U[G] . (1) 

Here iGQ 1 (x,y) = — (d 2 + m 2 )Sc(x — y) is the free Green's function and G is the full Green's 
function, both of which are defined on the closed time path C. The functional $[G] in ([TJ is 

1 However, their expressions are different from each other in the higher order terms of the skeleton expansion. 
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generally a sum of all possible 2PI graphs written in terms of G. A graph is called 2PI when it 
remains connected upon cutting two Green's function lines. 
The stationary condition for the effective action (fTJ 

gives rise to the Schwinger-Dyson equation for the Green's function G(x, y) 

G- 1 (x 1 y) = Gv 1 (x,y)-X(x,y) (3) 

with the proper self-energy defined as £ = iS$[G]/6G. The self-energy is divided into the local 
and the non-local part £ = £i oc + S non i- The £i oc contributes to the effective mass while the E n oni 
induces the mode-coupling between the different wavenumbers. The 2PI effective action should be 
invariant under the symmetry transformations of the system. Although we need to approximate 
the functional $[G] in practical applications, any truncation of $[G] which preserves the symmetry 
property gives the equations of motion consistent with the corresponding conservation laws |14[|15] . 

It is very useful to decompose the two-point function G(x, y) into two real functions, the sta- 
tistical function F(x,y) and the spectral function p(x,y) defined, respectively, as 

F(x,y) = \ <{0Or), #!/)}) = \ [G 21 (x,y) + G 12 (x,y)] (4) 

and 

p(x,y)=i([cP(x),cj>(y)}) =i[G 21 (x,y)-G 12 (x iy )] , (5) 

where (• • • ) represents the expectation value taken over a certain initial density matrix. The indices 
1 and 2 specify the branch of the contour C in the Schwinger-Keldysh formalism. The function F 
is called the statistical function because it turns out to be the Bose distribution function in the 
equilibrium state. The Schwinger-Dyson equation ([3]) can be equivalently rewritten in terms of 
Fix, y) and p(x, y) as coupled integro-differential equations 

(d 2 + rn 2 + j: ioc (x))F(x,y) = / dzT, F (x, z)p(z.y) - dzT, p (x, z)F(z.y) . (6) 

r° 

(d 2 +m 2 + Y, ioc (x))p(x,y) = dzE p (x,z)p(z,y) , (7) 

JyO 

where to is the initial time. Note that the non-local self-energy has been re-expressed similarly as 

Ms,!/) = \ [EwfolO+EJLfolO] , (8) 
E p (x,y) = ipi^^-S^^y)]. (9) 

The set of equations ((6]) and (0 is called the Kadanoff-Baym equation, which is the two-time 
formalism and describes the time evolution of the system from a certain initial configuration for F 
and p. Note that at each time step the spectral function p must satisfy the conditions following 
from the commutation relations: 

p(x,y)\ x o_, y0 = 0, 
d x op(x,y)\ x a^ y0 = (5 d (x-y), 
d x od y op(x,y)\ x a^ y o = 0. (10) 

Importantly, Eqs. ([6]) and (J7]) are non-local in time due to the so-called memory integrals 
appearing on the RHS. In other words, the evolution is non-Markovian depending on the evolution 
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Figure 1: Tadpole and sunset diagrams. 

history in the past. In many stable systems, however, the integrand of the memory integral dies 
away exponentially and the macroscopic time scale is separated from the microscopic one. 

It is instructive to consider the case of a uniform equilibrium state with a small value for the 
self-energy T, p (p°,p). Then we find that the spectral function p turns out to be the Breit-Wigner 
form (See Sec. [3]), 

P (Ap) = [(p0)2 _^_ sg/4 - ^<pW) a - , (ii) 

where il^ = p 2 + m 2 + RcS^ is the single particle energy including the mean- field effect. The 
arrow denotes the quasi-particle limit E p — ¥ 0. In this limit the p becomes a delta-function and for 
thermal equilibrium the statistical function F reduces to the Bose distribution 

F( P °,p) = 2n5(( P y ~ ^1) Q + ^Tl) • (12) 

In this paper we restrict ourselves to the spatially homogeneous situation. From the translational 
invariance, the statistical function F(x,y) — F(x a 1 y°,x — y) and the spectral function p{x,y) = 
p{x°, y°, x — y) can be Fourier transformed to F(x° ,y°;p) and p(x° ,y°;p). Then KB equations 
are simplified in the momentum space as 

(d 2 +p 2 +m 2 + S tad (x ))F(x , y°; p) - / dz°X F (x°, z°; p)p(z°, y°; p) 

r° 

dz X p (x ,z°;p)F(z ,y ;p), (13a) 

J to 

r° 

(d 2 + p 2 + m 2 + S tad (x°))p(x , y°, p)=- dz°E p (x°, z°; p)p(z°, y°; p) . (13b) 

Jy° 

Regarding the functional $[G], we approximate it with the skeleton diagrams obtained at the 
next-leading order in A. The self-energy E then becomes the sum of the local tadpole diagram and 
the nonlocal sunset diagram (Fig.[TJ: 

E loc (x) = E tBd (aO = ~F(x, x), (14) 
X«l(x,y) = ~G ab (x,y) 3 (15) 
where indices a, b denote the branch 1 and 2 of Schwinger-Keldysh contour C. Furthermore the 
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nonlocal part is divided into £ p and S p and written explicitly in terms of F and p as 

A 2 f d d k d d q 



Z F (x u ,z v ;p) = - — 



6 J (2ir) d (2tt) 



^/jp-k-q) 
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i^ , z°; k)F(z°, z°; q) - ^(z , z°; k)^ , z°; q) 



(16) 



S p (x ,z°;p) 



d d k d d q 
(2n) d (2n) d 



P(x°,z°;p k q) 



FOr , z°; k)F(x°, z°; q) - -p(x°, z°; k)p(x°, z°; q) 



(17) 



We solve these KB equations (fT3")) with the self-energy functions (fT4")l , (fT6"|) and (|17l) numerically in 
SeclU 



We need to specify the initial condition for p and F at x 



y 



t in order to solve this 



evolution equations. For the spectral function p it is fixed by the commutation relation as given 
in Eqs. (flU)) . For the statistical function F, we choose to set the initial conditions of the following 
functional form 



F(x°,y°;p) 
d x oF(x°,y°;p) 

.» „ EV„0 „.0.„\ 







x°=y°=t 


w(p) ( 


x°=y =t 


= o, 


x o =y o=t a 


= "(p)( 



1 



(18) 
(19) 
(20) 



where w(p) 2 = p 2 + m 2 and n p is a function we can freely specify. This form is assumed in analogy 
with the equilibrium solution in the quasi-particle limit. 

At later times in course of the evolution, we define the particle number distribution n p (X°) and 
the frequency w p (l°) in terms of F(x°, y°; p) [HI 123 Ell E3 Q] 



"p(^°)+2 



d x odyoF{x u ,y u ;p) 
d x od y oF(x ,y°;p) 



x°=y°=X t > 
1/2 



F(X u ,X u ;p)- \d x oF{x»,y u ;p) 



x =y =X 



o =v o =x o 



F(X°,X°;p) 



1/2 

(21) 
(22) 



Strictly speaking, these definitions (|2"Tj) and ([2"2"jl are valid only when the quasi-particle picture 
works well. Nevertheless, we expect that these quantities are good estimators to characterize the 
behavior of the system evolution. The system is expected to have a quasi-particle spectrum for a 
sufficiently small coupling A as shown in 1 + 1 [3TJ [35], 2 + 1 [23] and 3 + 1[2U [35] dimensions. 

Before proceeding to the next section let us compare the KB equations with the Boltzmann 
equation in 1+1 dimensions. In a homogeneous system the Boltzmann equation becomes 



d 



A 2 fd d Pl d d p 2 d d p 3 



1 



4 J (2ir) d (2ir) d (2ir) d 8n pi n p2 fi ps 
x K 1 + "psX 1 + n p )n pl n P2 - n P3 n p (l + npj(l + n P2 ))} 

x (2^) d+1 ,5 d (p 1 + P2 p 3 p)S(n pl + n P2 - n P3 - n p ) 



where fi p = \/p 2 + p?(t) and the mass p 2 (t) is the self-consistent solution of 

A f d d k n k (f) 



p(t) z = m z + 



(2n) d v / M (t)2 + k 2 



(23) 



(24) 
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In fact, this Boltzmann equation can be derived from the KB equations at the leading order in the 
gradient expansion and with the Markov and quasi-particle approximations 31 . We remark here 
that in 1+1 dimensions the Boltzmann equation cannot lead to thermalization because the particle 
momenta must be unchanged in each 2-to-2 collision in order to satisfy the energy and momentum 
conservations. 



3 Entropy of the relativistic Kadanoff-Baym equations 

In this section we derive the expression for the relativistic entropy in terms of the two-point functions 
G{x, y) for the theory, as an extension from the non-relativistic entropy current given in |27] 
and [28]. 

We start with the Schwinger-Dyson equation ([3]). Multiplying G from the right and left hand 
sides of Eq. ([3]), respectively, we obtain 

d 2 x + m 2 + ^G aa {x 7 x) 



dl + m? + ±G hh U,.!,) 



G ab (x, y)-ij dzXZA*, z)c cd G db (z, y) = ic ab 8{x - y) , (25) 
G ab (x, y)-i [ dzG{x, z) ac c cd Ef onl (z, y) = ic ab 5{x - y) , (26) 



where a and b assign the branch 1 and 2 of the Schwinger-Keldysh contour C and c ab = diag(l, — 1). 
We introduce the "center-of-mass" coordinate X = (x + y)/2 and the relative coordinate x — y. 
Then making the difference of these equations (|26|) and (|25j) and performing the Fourier transform 
with respect to the relative coordinate x — y, we find 



d i A f d d+1 k fdG aa (X,k) ( dG bb (X,k) 



d_ 

dp 



G 



ab 



dX 2 2 J (27r) d+1 V dX dX 
l jd{x- y )e*-C*-") / dz(E- nl (x, z)c cd G db (z, y) - G ac (x, z)c cd £f OJll (*, y)) , (27) 



where p and k are the momentum conjugate to x — y. When we make the sum of them and perform 
the Fourier transform, we get the expression 



' , 2 A ( f d d+1 k 
p — m — — 



J^yITl(G aa (X,p) + G bb (X,p)) 



G ab (X,p) 



ic ab + l - / d(x - y)e-^-^ / dz&ZJx, z)c cd G db (z, y) + G ac (x, z)c cd E^ nl (z, y)).(28) 



Starting the evolution at x° — y° — 0, we have G{x,y) only in a finite region of x° and y . It is 
therefore important to note that the interval of x° — y° is inevitably limited within ±X° in the 
Fourier transformation. 

The gradient expansion with respect to the center-of-mass coordinate X is adequate when the 
X-dependence of the system is smooth enough (See for example [33l [27]). We keep just the first 
order terms in the gradient expansion of Green's functions and the self energies here. For the 
expansion of the right hand side of Eqs. (|27p and (|28|) . we use the formula for two point functions 
K(x,y) and L(x,y): 



d{x - y)e lp - (x - y) / dzK(x,z)L(z,y) = K(X,p)L(X,p) 



dK dL dK dL 



2 \ dp^ dX,, dX^ 



O 



( d 2 



\dX< 



.(29) 



where K(X,p) and L(X,p) are the Fourier-transforms in x — y. We remark here the scale separation 
between X° and x° — y°. We implicitly assume that the time dependence on the former is smooth 
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and mild while the time correlation in the latter is much shorter. At the very early time of the 
evolution, this separation cannot be expected and the gradient expansion should be invalid. 

The derivation of the entropy current is most simplified in terms of the retarded propagator 
Gr = i(G n — G 12 ). To the first order in the gradient expansion, Eqs. (f2"T|) and (|2"5)) reduce to the 
equations for Gr: 



d(M — ±E P ) d d(M 



d 



G R (X,p) 



0. 



(M-^ p )G R (X,p) = -1, 



where M denoted 



M = p 2 -m 2 - ReS fl , E fl = E 11 - S 1 



In deriving Eqs. (|30f) and (|31j) we have used the well-known relations E 



11 



-22 



il2 



(30) 
(31) 

(32) 

V21 



G 1 



G 22 = G 1 



G and 2iImEj? = E„. The formal solution of the above simultaneous 



equations ([30]) and (f3Tj) is written as [33] 

G R (X,p) 



M - 



2^P 



(33) 



One should note here that M (E p ) is real (imaginary). Therefore, the real and imaginary parts of 
the retarded propagator Gr are given as 



ReG R (X,p) 
p(X,p)=2timG R (X,p) 



M 



M 2 ±E 2 



M 2 



IV2 



(34) 
(35) 



We see that the spectral function p(X, p) has the Breit-Wigner form (fTT]) in the first order approx- 
imation of the gradient expansion. 

Now we are ready for writing down the entropy current. The "derivation" goes somewhat in a 
heuristic way. We make the difference of Eq. ([27)) for (a, b) = (1,2) multiplied by ln(iG 12 /p) and 
Eq. (|2"T| for (a, b) = (2,1) multiplied by ln(iG 21 /p). Then we integrate the resultant expression 
over d d+1 p/ (2-K) d+1 to arrive at the following equation: 



9ns" = - 
2 



d d+i 



(27T) 



d+1 



£12 
G 21 



In— C(X,p). 



Here the term C, 



C(X,p) = i F(X,p) - Z F (X,p) p(X,p)) . 



(36) 



(37) 



may be identified as the collision term in the Boltzmann limit. With Eq. (|36|) . we define the entropy 
current s^(X) as 



d d+i 



P 



l^eG R 



V 



1 dReY, R 
2 



d 
dpt. 



dpp, 
E p iG 12 



In 



G 12 ln- 

iG 12 



P 
d 
dp^ 



G 21 ln- 



E p iG 



i p 

11 /^22\ 



2 ^ , iG 21 
In 



(38) 



where we have used the relations «(E n — E 22 ) = 2Re£^ and i(G 11 — G 22 ) = 2ReG_R. We have also 
applied the approximations E 12 ~ v 1 " V ' J v '' 
[27],|32]. 



Z„— and E 21 



Yip^y in the first order gradient expansion 



2 Please don't confuse this M with a mass function. Both M and £ have mass-dimension 2. 
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When we write the two-point functions in the form of the Kadanoff-Baym Ansatz G 12 = —ipf 
and G 21 = — ip(l + f) with a real function /, the above expression for the entropy current becomes 



^{f( P »-^)(-/W + ( 1 + /»M 1 + /»» 

4- ( -/Vn / + 4- ( "(I + f)) HI + f) 



dp, 



dpv 



After integration by parts over p M in the second line, we obtain a simple expression: 



d d+i 



P 



(2ir) d + 1 

where we introduced the notation 



P ~[p» 

i \ 2 dp ft 



ldReEfA T, p ldReG R 



i 2 dpn 



a(X,p) = -f\nf + (1 + f)ln(l + f) 



(39) 



(40) 



(41) 



One must distinguish this "occupation number" function / in the Kadanoff-Baym Ansatz from the 
distribution function n p defined in (12ip . 

Substituting the solution (|34[) for Gr, we can write the entropy current more explicitly as 




El 

4 



2M 



(2n) d + 1 

This expression further simplifies with use of ([331) tjl 




1 dRcSfl 

2 dp. 



P p dp^ 



4/ 



v 2 



a. (42) 



d d+1 p p 2 Y, f 
{2tt)< 1 + 1 2i 



P' 



2~dp~) ~ 2^% 



(43) 



This is one of the main results of this work. This expression of the entropy current is a natural 
extension to the relativistic case. The only difference between non-relativistic [37J [53] and our 
relativistic case is the factor | in front of the momentum derivative of the self-energy. We remark 
here that there is a discussion about the memory correction terms to the kinetic entropy in the 
non-relativistic case in Refs. [37] and [35] when we deal with the skeleton diagrams S non i beyond 
the NLO in A. 

In the quasi-particle limit, E non i — y 0, we know that G 12 = —ipf = 2nS((p ) 2 — flp )(#(— P°)+n p ) 
and G 21 = — ip(l + f) = 2ttS((p ) 2 — Cl 2 } )(8(p°) + n p ). In this limit the expression of the entropy 
current for p = reduces to the well-known form of the entropy density for bosons 



s° = 



d d p 
(2n) d 



[-n p In n p + (1 + rip) ln(l + n p )} , 



(44) 



as it should be. In general cases, however, the spectral function p(X,p) is defined as the Fourier 
transform of Eq. (JSJ) in x — y, and the occupation number function f(X,p) is then obtained with 
G 12 (X,p) — —ip(X,p)f(X,p). Although p°p( X, p) > in equilibrium, we have only a finite support 
in a; — y° in the initial value problem and the resultant Fourier transform p°p(X 1 p) may oscillate 
in p°, as shown in the next section. Accordingly the function f{X,p) can become negative, which 
brings a difficulty in evaluating the entropy density s° obtained at the leading order of the gradient 
expansion. 

Finally we show the fact that this entropy current obtained in the NLO in A formally satisfies 
the H-theorem. Namely, the RHS of Eq. (|36|) is positive semi-definite. This can be verified by 



The tadpole part should be the renormalized one in this expression in 1+1 dimensions. 
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Figure 2: Evolution of the distribution func- 
tion n p (o)p/m) from the tsunami initial con- 
dition (A/to 2 = 4). 



Figure 3: Evolution of the distribution func- 
tion n p (u]p/m) from the WS initial condition 
(A/to 2 = 4). 



substituting the expressions for £ p (fl~6|) and E p (Tl7|) into the RHS of (|36|) . As a result we obtain 
the relation 



d d+1 p 1 
(2 7 r) d + 1 2 n G 21 
1 A 2 /• d d+1 p d d+1 <j d d+1 A: d^ 1 ? 
8 ' 3T 



T (27r) d+1 (p + g- 



ln 



(27r) d + 1 (2Tr) d + 1 (2 7 r) d + 1 (27r) d + 
G 12 (p, X)G 12 (g, A)G 21 (k, A)G 21 (r, X) - G 21 (p, X)G 21 (g, X)G 12 (k, A)G 12 (r, X) 
G 12 (p, X)G 12 {q, X)G 21 (fc, X)G 21 (r, X) 



> 



(45) 



G 21 (p, X)G 21 {q, X)G 12 (k, X)G 12 (r, X) 

where we have used G 12 (— k, X) = G 21 (k,X). The last inequality holds since (x — y) In ^ > 0. Thus 
we proved that the H-theorem is fulfilled in the NLO in \<p A theory. However, at higher orders in 
the coupling constant A, the definition of the entropy current and the proof of the H-theorem are 
open problems. 



4 Numerical simulation 



In this section we show numerical results of time evolutions of the KB equations in the Atfi 4 theory 
in 1+1 dimensions. We discretize the space L = 2Na s into 2N grid points x n = na s (n 

X . —N+ 1, • • • , N — 1, N) with a s the lattice spacing and apply the periodic boundary condition. 
Accordingly the momentum has discrete value as p n = 2j ^ !l . The space derivative — d 2 is replaced 
with p 2 — -4- sin 2 , which removes much of the lattice artifacts [34] . We set N = 40, which 

is sufficient to study the momentum dependence. We also performed the simulation with N = 80 
and found no appreciable differences in the numerical results. We solve the evolution with the time 
step, a t /a s — 0.1. 

We set the mass ma s — 0.3 and varied the coupling A/m 2 =4, 2 and 1. We prepared the 
two different types of the initial conditions, "tsunami" distribution and the Woods-Saxon ( WS) 
distribution, respectively, 



< = ^ex P 



(H-Pt) s 
2a 2 



(46) 
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Figure 4: Evolution of the energy content in units of m from the tsunami initial condition with 
X/m 2 = 4; Kinetic e^i n /m (solid), tadpole e ta d/ m (dashed), sunset e sun /m (dotted) and the total 
energy e to t/"i (bold solid). See the text for the details. 



with a 2 /m 2 = 4.4 x (^) 2 , p T = 7 ■ 2ir/L and N T = 0.25, and 

n ws = _i 1 

P N W s e (v / IPl 2+m2 ~P»'s)/ K + 1 

with p ws /m = 2.936, n/m = 0.35 and M W3 — 0.5. The "tsunami" initial condition has two peaks 
at ±pt ("tsunami") and may be regarded as a toy model of the nuclear collisions. The WS initial 
condition is used to check the sensitivity of the evolution to the initial condition. These parameters 
in the WS case are tuned so that both the initial conditions give the same energy for X/m 2 = 4. 
Later in this section, we change the coupling constant A with other parameters fixed, in order to 
see how the evolution depends on the coupling strength. We monitor the energy conservation in 
each time step in order to estimate the numerical accuracy of our simulation. 

In Fig. [2] we show the time evolution of the number distribution n(Q p ) defined in Eqs. (|2"Tj) 
and (|22[) with the tsunami initial condition (|46[) . From this figure we confirm that our simulation 
reproduces the results of Ref. [21] . The peak of Gaussian distribution disappears rapidly and the 
values at high and low edge regions grow up with time. As a result, the particle number distribution 
approaches the Bose distribution function n p = l/(e ep//T — 1), with temperature T/m ~ 2.5 and zero 
chemical potential. Similarly in Fig. [3] we show the time evolution of number distribution n(Q p ) 
with the WS initial condition (|47|) . We see that n(ui p ) converges to the same thermal distribution 
as the one in the case of "tsunami" initial condition. 

Next we study the energy content of the system. The explicit expression of the energy is given 
in Appendix |A1 We plot the kinetic (|62a[) , tadpole (|62cp . sunset (I62d|) and total energy (|61a[) as 
a function of time mX° in Fig. [4J As for the total and kinetic energies we plot their differences 
measured from the initial value of the total energy e% t/m ~ 260 in our discretization. We find 
that the growth of the kinetic energy is canceled by the tadpole and the sunset energy to have a 
constant total energy. The energy is conserved within 0.5 percent in Fig. [4] 

4.1 Kinetic entropy 

Let us study the kinetic entropy (|4Q]) derived with the gradient expansion of the KB equation. To 
this end we first examine the shape of the spectral function p(X,p), which appears in the expression 
(|40j) and is needed to compute the occupation number function f(X,p) in a. The spectral function 
p(X,p) itself is physically important. In Figs. [5] and El we show p(X,p) for p x = 2nn/L with n = 0, 



(47) 
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Figure 5: Spectral function p(X,p°,p x ) with Figure 6: Spectral function p(X,p°,p x ) with 
p x = at mX° =1.5, 3 and 4.5. p x = 20tt/L at mX° = 1.5, • • • , 10.5. 
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Figure 7: Time evolution of kinetic entropy (|40p denoted in + and its quasi-particle approximation 
(|44p shown in a curve for the "tsunami" initial condition with coupling \/m? = 4. 



and 10, respectively, at several values of time X°, with the "tsunami" initial condition (|46|). We 
clearly see peak structures near p°/m ~ yj m 2 + p\ jra at later times in both figures. 

At early times, however, the spectral function p(X,p) shows oscillatory behavior. This can be 
understood as the uncertainty relation between the energy and the time. In the observation within 
a finite time interval — y°\ < X°, one can resolve thep dependence of p(X,p) on the scale larger 

than l/X°, because we have an oscillating factor due to J_ X( dt exp(— ip t) — 2 sin(p° X ) j p° . We 
numerically confirmed that the oscillation frequency is indeed proportional to X°. This means 
that any finer structure of p(X,p) than a scale 1/T is resolved only after the evolution time of 
X° > T. The sharper the peak structure is, the longer time it needs to be resolved, which is seen 
by comparing the cases with p x = and 20tt/L as shown in Figs. [5]and|6l respectively. 

The applicable range of the kinetic entropy [27l [28] has not yet been examined so far with using 
the numerical solution of the KB equation before this study. Given the oscillating p(X 7 p) near the 
initial time, we also have an oscillation for the occupation number function f(X,p) accordingly. 
We thus encounter a problem in evaluating a as it contains the logarithm of f{X,p). Here we 
come to recognize that the form of the kinetic entropy (j40[) obtained in the leading-order gradient 
expansion cannot be applied in the early stage of the initial value problem, although the gradient 
expansion itself becomes unjustifiable there. 

In Fig. [7] we try evaluating the kinetic entropy (|40l) as a function of time X°. Crude as it is, 
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Figure 8: The ratio between kinetic entropy 
(40) and its quasiparticle approximation (44) 
for the tsunami and WS initial conditions with 
A = 4. 
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Figure 9: Time evolution of total number den- 
sity for the tsunami and WS initial conditions 
with A = 4. 



we simply neglect the contributions from the region of p° where p(X,p) has negative values as 
an exploratory estimate. We obtain the entropy that decreases until mX° ~ 10. This peculiar 
behavior is presumably stemming from omitting the negative p(X,p) contribution, which seems 
likely to overestimate the entropy. In the later stage, say mX° > 20, the kinetic entropy increases 
monotonically as time proceeds, which is expected from the fact that the kinetic entropy (HQ)) 
satisfies the H-theorem. 

4.2 Entropy in quasi-particle approximation 

As an alternative estimate for entropy of the system, we use the the quasi-particle (QP) approxima- 
tion (|44[) with the number distribution n p (X°) defined in Eq. (|2ip. The quasi-particle approxima- 
tion (44) may be reasonable because the spectral function is nicely peaked [3T] near w 2 ~ m 2 + p 2 
as seen in Figs. [5] and [H Since n p (X°) is defined locally at time x° = y° without the Fourier 
transformation, we have no computational difficulty even at the very early stage, in contrast to the 
kinetic entropy (|40[) . We should remark, however, that n p (X°) is here obtained using the solution 
of the full KB equation via Eqs. fl^T} and 

Although the QP entropy (|44p yields somewhat a smaller value as compared with the kinetic 
entropy (|40p. the evolution profiles of these entropies are quite similar to each other, except at the 
early times. Smaller value for the QP entropy (|4"4"|) may be related to the fact that it neglects the 
finite width of the spectral distribution over p°. With the QP approximation we see in Fig.[7]that 
the entropy production is concentrated at early times mX° < 20 and slows down at later times 
mX° > 20, approaching an equilibrium value. 

In Fig. [5] we show the ratio between kinetic entropy (40) and its quasiparticle approximation 
(44) for tsunami and WS initial conditions. Until mX° ~ 10 the ratio decreases because of the 
numerical artifact, explained in the above. In the range of 20 < mX° < 80 the ratio becomes 
constant for both initial conditions. In the late time region mX° > 80 the ratio starts to increase. 
This tendency can be explained from the behavior of time evolution of number density. In Fig. [3] 
we show the time evolution of number density for tsunami and WS initial conditions with A = 4. 
At the late time mX° > 80 with A = 4 the number density decreases continually. In addition 
the distribution function shown in Figs. [2] and [3] has the small change for mX° > 80. Due to the 
decrease of number density and small change for n p the entropy with quasiparticle approximation 
(44) increases weakly compared with the kinetic entropy (40). 

We compare the evolutions of the QP entropy for three values of the coupling, A/m 2 =4, 2 and 1 
with the tsunami and WS initial conditions, respectively, in Figs. [TUJ and QT] In these figures, we see 
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Figure 10: Entropy density s° /m in the quasi- 
particle approximation for the "tsunami" ini- 
tial condition with coupling X/m 2 =4, 2 and 
1. The exponential fit with (f48l) is denoted by 
a thin line in each case. 
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Figure 11: Entropy density s°/m in the quasi- 
particle approximation for the WS initial con- 
dition with coupling X/m 2 =4, 2 and 1. The 
exponential fit with (|48h is denoted by a thin 
line in each case. 



A 


7o 


s max -4 7 


7o 


s max A 7 


4 
2 
1 


0.24 
0.084 
0.027 


2.54 0.16 0.030 
2.48 0.10 0.0031 
2.39 0.17 0.0024 


0.14 
0.036 
0.0085 


2.54 0.14 0.030 
2.44 0.13 0.0031 
2.39 0.19 0.0011 



Table 1: Slope parameter 70 near X c 
(right) initial conditions. 



and parameters in Eq. (|48j) for "tsunami" (left) and WS 



that the larger is the coupling the faster is the entropy produced and saturated to the equilibrium 
value. In order to quantify the approach to the equilibrium value, we fit the entropy evolution with 
a simple functional form 

s(X°) = W - Ae-^ mX °^ , (48) 

where s m ax, A and 7 are parameters. We chose to fit the evolution in the regions 100 < mX° < 150, 
300 < mX° < 600 and 600 < mX° < 900 for A=4, 2 and 1, because the approach to equilibrium 
is slower for smaller A. The resultant values for the parameters are summarized in Table. 1. The 
parameters 7 and s max take the same values independent of the initial conditions both for the 
coupling constants A =4 and 2. The A dependence of the parameter 7 seems non-trivial, contrary 
to the A 2 dependence naively inferred from (|4"5j) . For A=l, our fit seems still sensitive to the 
initial conditions. At the later stage 7 might coincide for both initial conditions, but in this range 
mX° > 900 we have energy errors more than 0.5%, so that we stop our simulation. 

4.3 Particle changing processes 

As seen in the behavior of n(u) p ) in Figs. [2] and [3l the system evolves toward the equilibrium 
state. The entropy in the QP approximation is also written in term of n(u) p ). Here we study which 
microscopic process contributes to the change of the distribution function n(uj p ) in course of the time 
evolution. Although the distribution n(w p ) in Figs. [2] and [3] is computed using Eqs. ([21]) and (122ft . 
it seems instructive to evaluate the time derivative dn p /dX° in the quasi-particle approximation, 
which is given in Eq. (|66[) in Appendix [Bj Within this approximation we can clearly separate 
out the contributions of 0-to-4, l-to-3, 2-to-2 and 3-to-l processes. In Figs. [T2l and [T3l shown are 
the contributions of each microscopic process on the RHS of dn p /dX° Eq. (|66|) at the momentum 
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Figure 12: Microscopic contributions to 
drip/dX at L ■ p = 7 ■ 2tt as functions of mX° 
in the case of "tsunami" initial condition with 
A/to 2 = 4. By integrating from mX° = 
to mX° = 35 each contribution to Sn p is ob- 
tained to be 0.51, 0.62, -3.17 and -0.24 for 0- 
to-4, l-to-3, 2-to-2 and 3-to-l processes. 



Figure 13: Microscopic contributions to 
dn p /dX° at L ■ p = 7 ■ 2ir as functions of time 
mX° in the case of the WS initial condition 
with A / m? = 4. By integrating from mX° = 
to mX° = 35 each contribution to 8n p is ob- 
tained to be 0.54, 0.75, -0.86 and -0.26 for 0- 
to-4, l-to-3, 2-to-2 and 3-to-l processes. 



p n — 2irn/L with n — 7 for "tsunami" and WS initial conditions, respectively. We set A/m 2 = 4. 
The mode with p is included as one of the three particles in the l-to-3 process while it is formed 
from three particles in the 3-to-l process. 

Even in the quasi-particle approximation, the number changing processes are possible because 
of the finite memory time as was studied in Ref. [37] • We see that at early times the number 
changing processes 1 -H> 3 contribute as well as2f>2 scattering processes. One should recall that 
dn p /dX° = in the Boltzmann limit in 1+1 dimensions; even the 2-to-2 process is possible only 
when it keeps the particle momenta unchanged. 

As is seen in Fig. [2 the momentum of n — 7 locates near the peak for the tsunami initial condi- 
tion, and n(u) p ) decreases rapidly toward the equilibrium value. In the quasi-particle approximation 
shown in Fig. |T2] the 2-to-2 process contributes largely to this decrease, although other number 
changing processes are also operative at early times, say mX° < 35. For the WS initial condition, 
the momentum of n — 7 locates near the shoulder position of the distribution and n(uj p ) decreases 
with time as seen in Fig. [3] In this case, the 2-to-2 and other processes seem equally contributing to 
dn(u} p )/dX Q . However, the 2-to-2 contribution is negative and rather non-oscillatory while others 
are fast oscillating. 

One should bear in mind the limitation of our discussion here. Namely, these observations are 
in the quasi-particle approximation, while we evaluated n(uj p ) using the full solution of the KB 
equation which includes both the effects of the finite memory time and the spectral distribution 
p(X,p). The latter is missed in the estimate with the quasi-particle approximation. 



5 Discussion 

We have studied the time evolution of the A(f> 4 theory in 1+1 dimensions in the framework of the 
KB equations with the statistical and spectral functions, F(X,p) and p(X,p), as basic ingredients. 
In this framework we can take into account two kinds of the "offshell" effects; one is the finite 
memory time effect and the other is the non-trivial form of the spectral functions p(X,p). 

We have seen that the particle number distribution n(cD p ), which is defined with the numerical 
solution of the KB equations, converges to the Bose distribution in the time evolution. We also 
showed that both the kinetic entropy (|40|) and the QP entropy ([44]) increase in course of the 
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evolution. In order to understand the mechanism of the entropy production, we analyzed the 
microscopic processes within the quasi-particle approximation and found that the 2-to-2 scattering 
processes as well as the particle number changing processes are operative in the early stage of the 
evolution. It is important to note here that in the Boltzmann description in 1+1 dimensions, the 
particle distribution function never evolves, dn(£j p ) / dX° — 0, as we mentioned previously. 

It is not clear, however, which effect is essential to the system thermalization, the finite memory 
effect or the spectral function p(X,p). It would be interesting to compare the finite memory time 
simulations with and without the quasi-particle approximation in 1+1 dimensions, in order to 
estimate the importance of the spectral structure p(X,p) in thermalization. We leave this for our 
future study. 

We have introduced formally the kinetic entropy (l40|) associated with the relativistic KB equa- 
tions, to the leading order in the gradient expansion. In general, the gradient expansion applies only 
when the X° dependence becomes gentle toward equilibrium, otherwise the higher order terms can 
give substantial effects. Furthermore, we encountered a difficulty in evaluating the kinetic entropy 
P0|) . the leading term in the gradient expansion, because of the fact that the spectral function 
p(X,p) shows an oscillation in p° with the frequency proportional to 1/X . This is understood 
as the uncertainty relation between the energy and the time. This yields negative values for the 
occupation number function f(X,p), and invalidates the expression for the kinetic entropy (|40[) 
at early times. This has not been recognized until our attempt for the numerical evaluation. At 
sufficiently later times when p(X,p) and f(X,p) become positive definite, the resultant entropy 
monotonically increases in time. In contrast, the QP entropy (1441) is expressed in terms of the num- 
ber distribution n p obtained from the two-point function F without the Fourier transformation, 
and it shows a monotonic increase in Fig. [7l Although the QP entropy has no strict foundation in 
the context of the KB equation, its time evolution seems physically quite suitable as an indicator 
of the system entropy. 

We have also analyzed and compared the behavior of kinetic entropy (|40[) and its quasiparticle 
approximation (]44|) . In the middle range of thermalization the ratio of both entropy becomes 
constant, so that entropy (|44p is a useful indicator of the thermalization. However at the asymptotic 
stage of thermalization due to the small change of n p entropy (|44|) is affected by the behavior of total 
number density. Therefore entropy (|44p has difficulty a little in estimating asymptotic behavior of 
the kinetic entropy (|4"4"|) which is based on KB equation, but does not take so large difference. 

6 Summary and outlook 

We have extended the kinetic entropy current to relativistic case from non-relativistic one on the 
basis of the KB equation \cf> 4 theory. The derivation was done in parallel to the non-relativistic 
case to the leading order of the gradient approximation and to the NLO of the skeleton expansion 
of $[G]. The derived kinetic entropy satisfies the H theorem. 

The numerical simulation of the KB equations in 1+1 dimensions has been performed with the 
"tsunami" and Woods-Saxon initial conditions. We have seen that the particle number distribution 
approaches the Bose-Einstein distribution with time, irrespective of the applied initial distributions, 
which reconfirms the results of Ref. [3T]. We noted that the Boltzmann equation with 2-to-2 
scatterings in 1+1 dimensions never thermalizes. Therefore, we see that the "offshell" effects 
included in the KB approach must be very important for thermalization problem. 

In evaluating the kinetic entropy (|40|) numerically, we recognized that the Fourier transforma- 
tion within the limited time interval for mX° < 10, makes the spectral function p(X,p) and the 
occupation number function f(X,p) oscillating in p°, accordingly the expression (|40[) must become 
complex- valued. At later times, on the other hand, the kinetic entropy (|40[) increases in time, indi- 
cating that the system approaches the equilibrium state. Therefore we find that the kinetic entropy 
(|40|) should be useful when the shape of the spectral function p(X,p) is well resolved, provided that 
the gradient expansion is also valid. 

As an alternative, we have studied the entropy evolution in the quasi-particle approximation 
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(|44[) . and found that it increases almost monotonically saturating at the equilibrium value. The 
late time behavior is similar to that of the kinetic entropy. We further studied the dependences of 
the evolution on the coupling constant and the initial condition. 

The KB equations involve the offshell effects of the finite memory time as well as of the non- 
trivial spectral functions p{X,p) and F{X,p). We examined effects of the finite memory time 
within the quasi-particle approximation, and showed that the 2-to-2 process as well as the number 
changing processes are operative for producing the entropy at early stage of the evolution. In order 
to study the importance of the spectral function p(X,p), it will be useful to compare the simulations 
with and without the quasi-particle approximation in 1+1 dimensions [3 7). We leave this for future 
study. 

The KB equation is one of the most promising approaches to describe the non-equilibrium pro- 
cesses in the quantum field theories. Although we limited our simulations in 1+1 dimensions to see 
the importance of the "offshell" effects, the simulations in high dimensions are certainly needed, 
where those effects will be likely important. In the field of high-energy heavy-ion collisions, early 
thermalization is one of the important issues under debate. The KB dynamics may provide a suit- 
able framework to investigate the early time evolution of such energetic nuclear collisions, and the 
entropy introduced in this paper may be one of the useful quantities to characterize thermalization 
of the system in course of its time evolution. Toward this end, extension and practical development 
to the case of gluodynamics are much desired. 
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A Energy Momentum Tensor 



We derive the expression for the energy-momentum tensor, <d^ u , of the \<j) A theory. The 2PI 



effective action T[G] Eq. (jTJ) is invariant under the translation 



x^ x H 



e M . Following Nocthcr's 



procedure, we apply the position dependent translation x^ — > x M + e M (x) to compute the change of 
the action <5r = J d v [e M (a;)0 A "' '(x)]. For e„ independent of x, we can prove the current conservation 
d v Q^ v — as a result of the invariance of the action ST = 0. The energy-momentum tensor 
reads from ST as the coefficient factor of d„e r (x) [15 | [29 j [24] . 

Under the position dependent translation, Green's function changes as 

G(x, y) -> G'(x, y) = G(x + e(x),y + e(y)) = G(x, y) + e x (x)dfG(x, y) + e A (j,)^G(x, y). (49) 

Then the change of each term in the action (JT]) is calculated as follows: the first term in Eq. ([T]) 
leaves no e term 



-TrlnG" 



G~ l {x, y) [e»(y)d]>G(y, x) + e"(x)d*G(y, x)] 



i / e"{x)d x J{x - x) = . 



(50) 



The second term gives rise to 



n—^-n 

~~ V .. ' / 1 ( ' 



[(dl+m 2 )5{x~y)} SG(x,y) 



>*{x)dl S(x - y)d*dl (G(x, y) + G(y, x)) 



-S(x - y)g lxv d£dlG(x ) y) + m 2 g^S(x - y)G(x, y) 



(51) 
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where we have used 



d£[5{x-y)G(x,y)] 



x.y 



5{x-y) [dZG(y,x) + dZG(x,y)] 



x.y 



The third term yields 
1 



5$ 



2j x<y 6G(x,y) 



(e»(x)dZG(x,y)+d»(y)dyG(x,y)). 



(52) 



(53) 



The third term 6 [^<E>[G]] can be rewritten more conveniently by use of the Jacobian. For example 
the tadpole part in $[G] changes under the translation as 



A / G{x,xf 



= 2A / G(x,x)(e»(x)dZG(x,y) + e(yrdyG(x,y)) 



(54) 



However we observe that the translation can be dealt with the change of variables in the integral 

as 



A / G(x + e(x), x + e(x)) 2 

J X 

and therefore the change can be recast to 



det 1 



dx» 
dx' v 



G(x', x 1 ) 



l\2 



A / G(x, xf 



-A 



, dx^ 



G{x',x') 2 = -A 



dxf* 



G(x, x) 



(55) 



(56) 



In this way, the total change of $ [G] can be re-expressed by using the Jacobian. Furthermore, the 
number of integration coincides with the power of the coupling A in general. This means that the 
Jacobians under the translation can be absorbed in the change of the coupling A — > \C( X ) an d that 
<5$ can be rewritten as fl5l [29l f24l 







[V 




2 





5$ 



5C(x) 



C=i 



(57) 



Since ST — J d u e^(x)Q tJ " 1 ' (x), the energy-momentum tensor is found as 



®^{x) 



5{x - y)d£8Z (G(x, y) + G(y, x)) - S(x - y)g^d^lG(x, y) 



+g> M "5(x-y)m 2 G(x,y) 



V— 

r sc(x) 



(58) 



In the case of uniform space, by taking the Fourier transform with respect to the spatial relative 
coordinate, we obtain the explicit expressions for the energy 00 and the pressure 8 11 afl: 



e 00 (X°) = \ /tJW +m 2 +d xO d y o)F(x , y °;p) 



15<f> 

— -^r— — counter term , (59) 

rO=yO =X n 2 S( ' v ' 



n (X^ 



d d p 



(27T) 



1 (p 2 -m 2 + d x od y o)F(x 7 y°;p) 



1 <5$ 



counter term . (60) 



Here the counter term cancels out the divergence in ^ . In 1 + 1 dimensions the divergence is only 
in the tadpole part which can be in the similar manner to that of [21] . 

In the NLO in A of the skeleton expansion, the tadpole and sunset diagrams contribute to 
8&/5(. Then the total energy and pressure are decomposed as 



Etot(X°) — E^i^X ) + i?tad(^°) + ^sun(^°) 
P tot {X°) = P kia {X°) + Ptad(*°) + Psun{X°) , 



(61a) 
(61b) 



The counter term in energy and pressure is the same. 
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where 



E kin (X°) = ~ J (0j(P 2 + ™ 2 + d x od y o)F(x°,y°; p)\ x o =y0=xo , (62a) 
P^n(X°) = l -J ^(P 2 -m 2 + d x od y0 )F(x°,y°; p)\ x o =y0=x o , (62b) 

£tad(*°) = -Ptad(^°) = \ J^- d ^ A {X°)F{X\X^p) + \ J^Srnl* F(X°,X°;p) , 

(62c) 



E san (X°) = -P sun (X°) = ~ |A^7(JC°,p) . 
Here Sm^ ad denotes the mass counter term: 



(62d) 



^ = -2/(2^2^' (63) 



j A f dd P 1 



and 



I(X°, p) = j x dt> [H p (X°, t'; p)F(t>, X°; p) - £ F (X°, f ; p)p(t', p)] (64) 
is finite in 1+1 dimensions @. 

B Microscopic processes in the Kadanoff-Baym equation 

The evolution of the particle number distribution can be understood in terms of the microscopic 
processes in the quasi-particle approximation. See [1] [37] for a similar discussion. 

First, we differentiate the distribution n p defined in Eq. (|2Tj) with respect to the time t to find 



- + rip ) d t n p = 



f dz°{ [E„(t, z°; p)F(z°, t; p) - E F (t, z°: p)p(z°, t; p)] $F(t, p)| t=t , 
•/to L 

- [£„(*, z°; p)d t F(z°, t; p) - £ F (i, z°; p)cV(A t; p)] F(t, f; p)}. (65) 



Meaning of the memory integrals on the RHS of Eq. (|6"5"|) becomes clear if we take the quasi-particle 
limit for p and F as given in Eqs. ([TTj) and (TT2"j) . We find t'; p) = w(p) _1 sin[(tD(p)(f — i')] and 
F(t, t'; p) = cl^p)" 1 cos{u)(p)(t — t')](n p + i). We substitute these to the self-energy E p and £ F at 
the NLO in Eq. (lB5"j) and arrive at the following expression: 

„ f \ A 2 f d d q_ d d k f\. 1 

dtn P (t) = - / 7^7^-, / dt> 



3 7 (2?r) d (2ir) d J to 2J)(p)2w(q)2J)(k)2w(p - k - q) 
| [(1 + n p )(l + n q )(l + n k )(l + n p _k-q) - ri p n q n k np_ k - q (0] 
x cos [(cD(p) + cD(q) + w(k) + cD(p - k - q)) (i - i')] 

+3 [(1 -I- n p )(l + nq)(l + n k )n p - k - q - n p n q n k (l + n p _ k -q) (*')] 
x cos [(w(p) + w(q) + w(k) - cD(p - k - q)) (i - i')] 

+3 [(1 + 7ip)(l + rt q )n k n p _ k _ q - n p riq(l + n k )(l + ra p _ k _ q )(i')] 
x cos [(w(p) + w(q) - w(k) - Q(p - k - q)) (t - t')] 

+ [(1 + ^ P )^q^ k n p - k -q - n p (l + n q )(l + n k )(l + n p _ k _ q )(t')] 
x cos[(d)(p) -w(q) -w(k) - w(p-k-q)) (t-i')] }• ( 66 ) 



'In the case of 3 + 1 dimensions, the renormalization of self-consistent theories is required and developed in 
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We can interpret physically each microscopic process in the bracket {• • • } of Eq. (|66|) . The first term 
represents creation and annihilation processes of four particles 0f)4. The second term describes 
the process 1 — > 3 where n p is involved as one of the three particles and its reverse. The third term 
corresponds to 2 ■<-» 2 scattering process. The last term describes the decay of n p to 3 particles and 
its reverse. 

The number changing processes are allowed because of the finite memory time t — t'; the energy 
conservation in each microscopic process can be violated [37] . We evaluated these contributions and 
showed then in Figs. [T2l and [T51 Oscillatory behaviors seen in these figures obviously come from cos 
functions. Removing the memory time effect by taking the infinite time limit t— t — > oo, we recover 
the strict energy conservation in the microscopic process, \im t -t ->oo J t dt' cos(uj(t — t')) = ■k5{lo). 
In this limit only the 2 2 scatterings are allowed and then we obtain the Boltzmann equation 

(23). 
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